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We show how expansions in powers of Planck's constant h — h/2-n can give new insights into 
perturbative and nonperturbative properties of quantum field theories. Since fi is a fundamental 
parameter, exact Lorentz invariance and gauge invariance are maintained at each order of the expan- 
sion. The physics of the h expansion depends on the scheme; i.e., different expansions are obtained 
depending on which quantities (momenta, couplings and masses) are assumed to be independent 
of h. We show that if the coupling and mass parameters appearing in the Lagrangian density are 
taken to be independent of fi, then each loop in perturbation theory brings a factor of h. In the 
case of quantum electrodynamics, this scheme implies that the classical charge e, as well as the fine 
structure constant are linear in h. The connection between the number of loops and factors of h is 
more subtle for bound states since the binding energies and bound-state momenta themselves scale 
with h. The h expansion allows one to identify equal-time relativistic bound states in QED and 
QCD which are of lowest order in h and transform dynamically under Lorentz boosts. The possibil- 
ity to use retarded propagators at the Born level gives valence-like wave-functions which implicitly 
describe the sea constituents of the bound states normally present in its Fock state representation. 

PACS numbers: ll.15.Bt, 12.20. Ds 



I. INTRODUCTION 



Planck's constant % = h/2'K is the fundamental constant of nature related to quantum effects [1]. The most familiar 
application of the Planck constant is the commutation relation [xi,pj] = ihSij vifhich limits simultaneous measurements 
of position and momentum. More generally, h enters explicitly in the commutation relations of conjugate operators 
and fields, thus providing the fundamental basis of quantum field theory. Each order of an expansion in h must obey 
all the underlying symmetries of the theory. 

The Planck constant has units of action; i.e., the product of energy and length in units where the velocity of light 
c = 1. One commonly assumes units such that h = 1. However, it is illuminating to keep the occurrence of powers of 
h explicit since this allows one to distinguish quantum versus classical physics. For example, the AdS/CFT prediction 
for the ratio of shear viscosity to the entropy density of a multi-particle system has the lower limit [2] 77/s > H/Att. 
The origin of h in this relation can be traced to the assumed quantum form of the entropy of a black hole in a higher 
dimensional theory. The empirical observation that rj/s in heavy ion collisions at the Relativistic Heavy Ion Collider 
[3] is not far from the AdS/CFT prediction thus suggests that the dynamics of high energy central heavy- ion collisions 
is in the quantum domain. 

Physical phenomena are usually expected to follow the laws of classical theory in the (hypothetical) limit h ^ 0. 
Surprisingly, this is true only for a careful choice of /i-independent quantities (momenta, couplings and masses). A 
simple illustration is provided by the standard harmonic oscillator in nonrelativistic quantum mechanics where the 



potential is V{x) 



The propagation of a particle from {ti,Xi) to {tf,Xf ) is given by the path integral 



ti) = j [Vx{t)] exp 



2h 



dt{x^ - io^x^) 



= / exp 



~2 



(1) 



In the second equality we have removed the explicit dependence on h by scaling the coordinates as ^ = xj^/h. (The 
scaling of the Jacobian is irrelevant for this discussion.) Remarkably, the full quantum mechanical structure of the 
harmonic oscillator model persists as — >■ when one uses the scaled variables ^, since the propagation — >■ ^/ is 
independent of fi, as are the scaled bound-state energies e„ = En/h = w(n -|- |). Thus there is a domain of positions 
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X oc and momenta mx oc \/h where the action 5* is proportional to h and the system stays quantum mechanical 
even in the /i — limit. On the other hand, the propagation between fixed (fi- independent) positions Xi,Xf involves 
large values of ^i,^/ oc 1/VH, and thus the transitions between highly excited levels (with n of order 1/h) correspond 
to classical dynamics in the h limit. 

In the case of general relativity, H can be eliminated [4] from the equations of motion, and physical phenomena 
only depend on dimensionless quantities such as uqed- Although mass cancels out of the equations of motion in 
classical gravity (due to the equivalence principle) , it appears in the Schrodinger equation [-5] , but always in the form 
m = m/h. As we shall show, to is the fundamental mass parameter which appears in the equations of motion for 
fields when one formulates quantum field theory through the action principle and functional integrals. 

In quantum field theory the T-matrix elements of lowest order in h are usually considered to be tree diagrams 
(Born approximation), with each loop correction introducing one additional power of h [6]. However, Donoghue and 
Holstein and their collaborators [7] have demonstrated that classical physics (of lowest order in h) can also emerge 
from loop diagrams where zero mass quanta appear. As we shall see, the difference arises from the definition of the 
S — )■ limit, i.e., the limit depends on which Lagrangian parameters are taken to be independent h. 

Born diagrams usually provide a good first approximation to scattering amplitudes in quantum field theory How- 
ever, the bound-state poles of a scattering amplitude are not present in tree (or any finite number of loop) diagrams, 
but instead are generated by the divergent perturbative expansion of the covariant Green's function. For example, the 
Schrodinger and Dirac bound states, which arise from tree-level interactions of an electron in an external Coulomb 
potential, emerge in field theory from the infinite sum of ladder and crossed-ladder Feynman diagram contributions 
to the electron- muon Green's function in the limit where the muon mass is taken to infinity [8] . The binding is caused 
by loop momenta which are oc h, thus changing the relation between the number of loops and the power of h. 

We shall argue that the h expansion can provide a systematic approximation scheme for bound states formed by 
interactions between particles. This expansion is equally valid for relativistic and non-relativistic dynamics. We 
find that the lowest order interaction kernel in h indeed defines a viable "Born term for bound states". This Born 
approximation is insensitive to the ie prescription of propagators, which allows a simple Hamiltonian equal-time 
development in cases where the Coulomb interaction dominates. A hidden Lorentz boost covariance provides a 
non-trivial test that the approximation correctly includes all contributions of lowest order in h. 

II. THE h EXPANSION 

There is a general understanding that each loop contribution to quantum field theory amplitudes is associated with 
one factor of h [6]. Loop integrals represent quantum fluctuations and classical physics is expected to emerge from 
the Born term in the /i — > limit. To fully define this limit one needs to specify the h dependence of the fields, 
couplings and mass parameters of the action S. The choice is not as obvious as it sounds since h appears not only 
in the path integral measure exp{iS/h) but (for dimensional reasons) also in the action itself. In this section we 
establish the scaling with h which makes the loop and h expansions equivalent. In specific problems the energies and 
momenta may also scale with H, as illustrated by the harmonic oscillator above. In the next section we show how the 
contribution from momenta of O (h) in the multiple loops of ladder diagrams describes the physics of scattering from 
a fixed potential and gives rise to Schrodinger and Dirac bound states. 

We use units where the speed of light c — eo — 1. Planck's constant h then has dimension of energy E times 
length L; i.e, [h] = E ■ L. The dimensions of the fields and other parameters of the action are as usual fixed by the 
requirement that S/h be dimensionless. 

The gluon part of the QCD action, 

-^J d'^xiG"")^ ^ - ^ J d^xid^'A" -^•'A^' + i^[A^',A•']y (2) 

sets [A] = E^/"^ ■ L^^l"^. The gluon coupling g has dimension [g] ~ E^/"^ ■ L^/^, ensuring that as — g-^/inh is 
dimensionless. This requires that g appear with an inverse factor of h in (2). For the same reason the classical electric 
charge e and mass to are divided by h in the action Ssqed of scalar QED, 

SsQED ^ j d^x{ [{d^ + I '-A^)<t>] ^ [{d^ + z ^-A^)<t>] - ^0t^) (3) 

The boson field dimensions [</>] = [A] = E^'"^ ■ L'^/^ implied by (2) and (3) agree with the dimensions of the 
corresponding classical fields. Fermion field dimensions are convention dependent since Grassmann numbers have no 
classical counterparts. While factors of H are rarely displayed in quantum field theory, the convention of Ref. [9] is to 
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write the electron part of the QED action in the form 

Sqed^ j d^x[i^{ih$-eA-m)ilj] (4) 

with a factor h in the kinetic term. This gives the electron field the dimension [tj}] = L"^^"^ so that ■0^''/' has the 
dimension of a probability density. With this convention the classical charge and mass appear in the action without 
a factor of /i, and the anti-commutator of the electron field is 

{^\t,x),^{t,y))=8\x-v) (5) 

Insofar as Ti is regarded as a constant of nature, the factors of h in the action only serve to ensure the proper units of 
the parameters. However, to define an % expansion we need to specify the dependence on h of all quantities appearing 
in the action. In analogy to the harmonic oscillator (1) we take the boson fields to scale as 

A = Aj^Jh, EE (6) 

with A and independent of h. Thus \A\ ~ [0] — L^^. Similarly, we take the parameters 

g=-, e=- and m=- (7) 

to be independent of h. Thus m has the dimensions of wavenumber [to] = L^^. In the functional integral measures 
for the gluon and scalar fields, 



^ SsQED = I d^x [{d^^ + I ~e^A^^)$\ ^ [{d^ + i ~e^A^)$\ - m^4>^ (8) 



H then appears exclusively in the combinations gV^ and e^/h. Hence loop corrections of O (g^) and O (e^) will be of 
Oih). 

With the convention of a factor h in the derivative as in (4) we take the fermion field to be independent of h, while 
the coupling and mass should scale as in (7). Then 



^ Sqed — Id''' 



$ — eVhj^ — m)%l! (9) 



again ensuring that each loop is of O [K) . 

The present derivation is equivalent to (albeit perhaps more transparent than) the standard one of Ref . [6] , which 
associates a factor h with each propagator, an hr^ with each vertex and assumes the parameters appearing in the 
action to be independent of h, i.e., g,e,rh (x hP as in (7). Thus the free gluon propagator (AA) = h{AA) (x h since 
{AA) oc at Born level. Similarly the Born term of (AAA) = h^/'^ {AAA) (x gh?. The fermion propagator {ipi') is of 
O (^°) with the field normalization implied by (4). Each loop correction adds one power of h.^ 

According to (7) the equivalence of the loop and H expansions thus, remarkably, requires that the classical couplings 
g, e and masses to divided by h are fixed in the /i — > limit. For example, the dimensionless fine structure constant, 

"=4^=4^ 

is then cx h. The alternative of keeping e and m independent of h was considered in Ref. [7], where it was found that 
classical physics can emerge from loops in an amplitude where pairs of massless photons appear in the i-channel. 

Even though variations of h can not be studied experimentally, an H expansion is useful as a formal tool which 
clarifies the structure of theory. The most appropriate expansion scheme depends on the physics at hand. 

In the remaining part of this paper we shall be mostly concerned with bound-state equations which require summing 
an infinite number of perturbative diagrams. The h expansion will allow us to identify a Born term for bound states 
to which loop corrections may be systematically applied. In particular, we shall make use of the fact that Born terms 
are insensitive to the ie prescription of the propagators, whereas that prescription is crucial for loops. 



Reference [6] uses the normalization convention for fermion fields where there is no factor of h in the kinetic part of the action (4). 
Then the fermion fields scale similarly to the boson fields (6), and identical results are obtained for the powers of h in amplitudes. The 
normalization convention only affects the Jacobian of the functional integral (and the external lines of the Green function). 
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FIG. 1: The first two ladder diagrams contributing to non-relativistic atoms in the hmit of smaU couphng a are shown. 



III. APPLICATION OF THE h EXPANSION TO BOUND STATES 



Bound states in quantum field theory can be identified from the divergence of the perturbative series of the T-matrix 
at each bound-state energy. In the familiar case of non-relativistic QED atoms, the divergence arises from ladder 
diagrams with loop momenta which scale with a. The binding energy of hydrogen thus scales as En/rrie — —■^a^/n'^ oc 
It?. Hence an h expansion at bound-state poles requires that the external momenta in the T-matrix scale with h. In 
such a limit the usual relation between the number of loops and power of h need not hold. The ladder diagrams are 
in fact sensitive to infrared loop momenta that contribute inverse powers of h. The lowest-order contribution in a 
to the hydrogen atom is also of lowest order in h, i.e., it is a Born term for the bound state. The loop momenta do 
not scale with h in the higher-order corrections to propagators and vertices, implying the usual relation between the 
number of such loops and the power of h. 

The bound-state poles of non-relativistic QED atoms are obtained by summing ladder diagrams, the first two of 
which are shown in Fig. 1. The tree diagram in (a) generally dominates the loop diagram (b) at small coupling strength 
a. However, for external momenta within the range of the atomic wave function the loop diagram is unsuppressed 
even in the a — > limit. In the center-of-mass system this requires \pi\ — \p2\ to be of O (am) {m being the mass 
of the constituents) and |p° — to| to be of O (a^m). The restricted range of the loop momentum, d'^k oc in the 
bound-state domain, as well as the additional factor of a in the loop diagram, is balanced (relative to the tree diagram) 
by its three additional propagators, each of which is of O (l/a^). The situation repeats for all of the higher-order, 
multi-loop ladder diagrams, allowing the ladder series to diverge and bound-state poles to occur at any value of a. 

Since a = e'^h/Air cx h the previous argument demonstrates that the loop diagram in Fig. 1(b) is, in the bound-state 
domain, of the same order also in h as the tree diagram of Fig. 1(a). The loop does not bring an extra factor of h since 
we consider a limit where the (external and internal) momenta depend on h, similarly to the case of the harmonic 
oscillator (1). Repeating the argument for ladder diagrams with more rungs, we may conclude that the Schrodingcr 
equation defines the Born term for non-relativistic atoms. The standard higher order corrections in a to the binding 
energies and wave functions are also of higher order in h. 

The energy fc" cx a^m exchanged in the ladder loops is of higher order in h than the 3-momentum |fc| cx am. Thus 
in the Born approximation we may set A:" = in the photon propagators. This makes the non-relativistic bound-state 
dynamics equivalent to tree-level scattering from an external, instantaneous potential. 

The rclativistic Dirac-Coulomb electron bound-states are obtained by summing e/i scattering diagrams which include 
all ladder and crossed- ladder photon exchanges [8] . In the limit where the muon mass is large, the sum of crossed and 
uncrossed photon exchanges gives a S{k'^) which suppresses energy exchange and reduces the dynamics to tree-level 
scattering. In this sense the rclativistic Dirac-Coulomb bound states are also of lowest order in h. 



IV. BOUND STATES IN A STATIC EXTERNAL POTENTIAL 



The ie prescription in propagators is related to the boundary condition in time and is irrelevant at lowest order in h. 
Heuristically, this is seen from the damping factors exp(±t5) which give energy denominators Einuiai— Eintermediate^i^ 
where e = 5h. In perturbative expansions of the S-matrix the tree diagrams are in fact independent of ie. Conversely, 
the ie prescription determines the discontinuities of diagrams which are of higher order in h through the pinching 
of loop integrals. For Born level bound states this prescription independence allows alternative but equivalent wave 
functions as we shall next discuss. 

Consider first the Dirac-Coulomb bound states in a static external potential. In four-momentum space the in- 
teractions with the Coulomb potential A^{k) do not change the energy component of the particle's momentum. 
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Denoting a single Coulomb photon interaction by K the Green's function G of the particle can be expanded as 

G(/,p) + SKS + SKSKS +... = S + SKG = ^i^'^} + ... (11) 

p-' — bj 

In the last equality we displayed the pole contribution of a bound state with energy whose residue R(E,p) is easily 
seen to satisfy the Dirac equation. Since K preserves p'^ the Green's function G is independent of the ie prescription 
at the negative-energy pole of the Dirac propagator S{p). In other words, for > we have p'^ + Ep > 0, where 
Ep = \/ + m? . In particular, the bound-state energy E is the same whether we use a Feynman Sp or retarded Sn 
propagator, 

„ f + m 

If we time-order the interactions through a Fourier transform p^ ^ t the bound state gives a stationary contribution 
G{t,p) = exp{—iEt)R{E,p) + . . . However, the Fock state decomposition of its equal-time wave function depends on 
the choice of propagator (12). In the Feynman propagator Spit-iP) the negative energy components move backward in 
time. This gives rise to "Z" -diagrams describing particle-antiparticle pair fluctuations in the time-ordered interactions 
with the static potential. The equal-time wave function of a Dirac-Coulomb bound state thus has Fock components 
with any number of pairs, and its explicit expression is, to the best of our knowledge, not known even in simple cases 
such as a l/r potential. 

In the case of the retarded propagator Sj^ (12), the negative-energy components move forward in time, 

SR{t,p) = ^ [(i?p7° -P-1 + m)e-^^-* + {Ep^"+p- J - m)e^^^'] (13) 

The corresponding time-ordered interactions have no Z-contributions, thus only the single bound particle is present 
at any intermediate time. As a consequence of this, the retarded propagator is local in space at infinitesimal times, 

lim SR(t.x) = 7° 5^(3;) (14) 

By contrast, the right-hand side would be non-local in x for the Feynman propagator Sp- This may be understood 
by invoking a completeness sum over states at an earlier time: the particle first moves backward in time and then 
returns to a;' 7^ a; at t = 0. 

The absence of pair production in retarded propagation as well as the locality property (14) allows a one-particle 
Hamiltonian description as in relativistic quantum mechanics [10]. The bound-state wave functions are then given 
by the Dirac equation. These single-particle wave functions describe the same bound states, which using Feynman 
propagators, contain an indefinite number of particle pairs arising from Z-diagrams. As we have seen from (11), the 
bound-state energies E of both pictures agree at lowest order in h. 



V. RELATIVISTIC BOUND STATES IN FIELD THEORY 



The Hamiltonian description of relativistic bound states in an external Coulomb potential can in certain cases be 
extended to field theory [10]. According to the above arguments the binding energies will at Born level be independent 
of the ie prescription used in the bound particle propagators. This allows the use of the retarded propagator Sr (13) 
which suppresses Z-diagrams and makes the time evolution local as in (14). Coulomb exchange gives an instantaneous 
potential due to the absence of OqA^ terms in the Lagrangian density. Transversely polarized photons propagate in 
time and thus they generate higher Fock states. Here we only consider cases where Coulomb exchange is dominant. 

In Z? = 1 + 1 dimensions only Coulomb exchange contributes in a gauge where — 0. For an e^//^ state where 
the electron is at xi and the muon at X2 

\E,t = 0)= J dxidx2^l{t = 0,xi)x{xi,X2)i^^{t = 0,X2)\0) (15) 

the equation of motion (Gauss' law), 

— dlA°{x; xi,X2) = e [6{x — xi) — S{x — X2)] (16) 
determines the instantaneous Coulomb field. 



- A° (x; xi, X2) = le[\x - xi \ - \x - X2\] 



(17) 
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It is then straightforward [10] to determine the condition on the bound-state wave function x{xi,X2) which ensures 
a stationary time development, 

l°{~idxt 7^ + ™e)x(a;i,X2) - x{xi,X2)l^{id 7^ + m^) = [E - V{xi - X2)\x{xi,X2) (18) 

Here the kinetic term of the electron operates on x(xi,X2) from the left, and that of the muon from the right. The 
potential following from (17) is 

V{x) = \e^\x\ (19) 

In _D = 1 + 1 dimensions the Dirac matrices as well as the wave function xi^ii ^2) niay be taken to be 2 x 2 matrices. 

Not surprisingly, the bound-state equation (18) has a "double Dirac" form, and as such was proposed already by 
Breit [11]. Now this equation is seen to provide an approximation of lowest order in h to relativistic bound states, 
with the potential (19) in l-fl dimensions fixed by QED. 

A stringent test that (18) actually represents the exact Born level result is that it is consistent with the Poincare 
invariance of QED. Translational invariance is explicit, but the boost invariance is hidden (dynamic), since the 
constituents are at equal time in all frames. It turns out [12] that for the linear potential (19) (and for no other form 
of the potential) the bound-state energy indeed has the correct dependence on the CM momentum P of the bound 
state, 

E = V'P2TPM2 (20) 

The P-dependence of the wave function xi^ii ^2) is non-trivial. It resembles an ordinary Lorentz transformation, but 
with E ^ E — V{x) (the canonical energy) in the boost parameter. Thus the wave function contracts at a rate which 
depends on the separation x between the constituents. 

The properties of the relativistic bound states defined by (18) merit further study. At large separations between 
the constituents, where V{x) 3> E^ the wave function has a constant, non-vanishing density in x. This may reflect 
the particle pairs which are polarized from the vacuum in a formulation using Feynman propagators. 

The two very different but equivalent pictures of bound states that we find here using retarded vs. Feynman time 
development may be related to the well-known puzzle of the quark model vs. parton model views of real hadrons: 
Hadron excitation spectra reflect mainly their valence quark degrees of freedom, even though sea quarks contribute 
importantly to deep inelastic scattering at low xbj and low [13]. 

A bound-state equation analogous to (18) may be derived in _D = 3 + 1 dimensions by assuming a non-trivial 
boundary condition in Gauss' law (16). The homogenous solution for then gives rise to a linear potential. Poincare 
invariance is respected similarly to the 1+1-dimensional case. It is also possible to derive the analogous Born level 
meson and baryon bound states in QCD [10]. The qqq potential is gauge covariant and confines the three quarks 
in a symmetric way. In some respects this resembles the soft- wall AdS/QCD models [14, 15] which utilize a linear 
potential in an effective Dirac equation in AdS space. 

VI. CONCLUSIONS 

The Planck constant h — hjl-K is a fundamental constant of nature which gives a measure of quantum effects. In 
this paper we have shown how expansions in powers of Planck's constant h can give new insights into perturbative and 
nonperturbative properties of quantum field theories. Since fi. is a fundamental parameter, exact Lorentz invariance 
and gauge invariance are maintained at each order of the expansion. 

It is common to set = 1 since there is a general belief that the power of h is given by the number of loops in 
the perturbative expansion [6]. In the functional integral h appears in the measure exp(i5/S), and the action S is 
usually assumed to be independent of h. It is then argued that classical physics emerges in the /i — >■ limit since the 
rapidly varying phase S jh selects field configurations for which the action S is stationary. This argument is, however, 
somewhat oversimplified since there may be field configurations which make S on h. We have illustrated this with the 
harmonic oscillator, whose full quantum mechanical bound-state spectrum persists in the /i — > limit. 

In fact, the physics of the h expansion depends on the scheme; i.e., different expansions are obtained depending on 
which quantities (momenta, couplings and masses) are assumed to be independent of h. We have shown that if the 
coupling and mass parameters appearing in the Lagrangian density are taken to be independent of h, then each loop 
in perturbation theory brings a factor of h. In the case of quantum electrodynamics, this scheme implies that the 
classical charge e, as well as the fine structure constant are linear in h. 

We have also noted that for dimensional reasons the QCD coupling g is divided by h in the action (2). Similarly, 
the classical charge e and mass m are divided by h in the QED action (3). 
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It may seem natural to fix the classical quantities e and m as fi, — >■ 0, as is done in [7]. However, this scheme 
introduces a non-trivial dependence on h into the action S. We also note that the fine structure constant a = /Anhc 
diverges in the S — >■ limit if the classical coupling e is held constant. Furthermore, it has been demonstrated through 
explicit examples in this scheme that loops containing massless quanta, such as the gravitational form factors of the 
electron, can contribute to classical physics in the h ^ limit [7]. The underlying reason is that an h expansion 
is not uniquely defined without specifying which quantities (momenta, couplings and masses) are to be regarded as 
independent of h. 

On the other hand, if the parameters of the action, g — g/h, e — e/h and rh = m/h are taken to be independent 
of h, the fine structure constant a — e^h/Air cx h and the h expansion is equivalent to the perturbative expansion. 
For example, in the fi. — > limit, the loops which define the /3 function vanish, so that the running coupling e(/^^) is 
constant as in a conformal theory. 

The h expansion is particularly non-trivial and illuminating in the case of bound states. We have considered whether 
one can define a Born term for bound states, which is equivalent to a tree (no loop) approximation when g, e and rh 
are fixed. This would allow an unambiguous and physically motivated approximation scheme for relativistic bound 
states, maintaining all symmetries of the theory at each order in h. In the familiar case of non-relativistic QED atoms 
the binding energy Eund oc a^m cx h^ih depends on h. In order to stay on the bound state pole of a Green's function 
in the /i — !■ limit we must allow the momenta to scale with h, which introduces additional sources of h. This is why 
the sum of multi-loop ladder diagrams in perturbation theory reduces to a Born approximation for atomic states. 

Born terms are insensitive to the ie prescription of perturbative propagators. This is explicitly seen for tree diagrams 
of the S-matrix, and we have shown that it holds also for bound states in a Coulomb potential. Using Feynman or 
retarded propagators does, however, make a consequential difference for the wave functions of the equal-time bound 
states. Relativistic bound states given by the Dirac equation have an indefinite number of particle pair constituents 
in their Fock expansion which arise from Z-diagrams due to Feynman propagators. In contrast, if one chooses 
retarded propagators, the bound state appears to contain just a single particle (of positive or negative energy) whose 
distribution is given by the standard Dirac wave function. 

Our understanding of bound states formed by the mutual interactions of relativistic particles is still very limited. 
This may be due to the lack of a physically well motivated and manageable first approximation akin to the tree 
diagrams of scattering amplitudes. In this paper we have shown that there is a well-defined Born approximation for 
relativistic bound states. With retarded propagation one obtains simple and specific bound-state equations. The 
two-body equations in QED and QCD have a hidden boost invariance which strongly suggests that they include all 
effects of lowest order in h. Just as for the Dirac equation, the resulting valence wave functions implicitly describe 
the multiple pair constituents generated by Feynman propagation. 
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